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Abstract 

This paper discusses the dynamical properties of p-spin models 
with Kac kind interactions. For large but finite interaction range R 
one finds two different time scales for relaxation. A first relaxation 
roughly independent of R where the system is confined to limited 
regions of the configuration space and an R dependent time scale 
where the system is able to escape the confining regions. 1 will 
argue that the R independent time scales can be described through 
dynamical mean field theory, while non-perturbative new techniques 
have to be used to deal with the R dependent scales. 



1 Introduction 

Mean-Field spin glass model, describe glassy phenomena as sharp ergod- 
icity breaking transitions jX,. It is well known that in that context two 
basic route to glassiness are found. In the SK and similar models, one 
finds second order phase transitions with a diverging susceptibility and a 
continuous order parameter I^. In p-spin models, one finds a transition, of 
the second order kind from a thermodynamic point of view displays a first 
order jump in the order parameter j2j. While the SK model was proposed 
as possible starting points to understand the physics of spin glass materials, 
p-spin like models have deserved a lot of attention as mean-field schemati- 
zation to understand the structural glass transition of supercooled liquids 
0]. In particular, the physical analysis of the model shows the presence of a 
phase where despite no sign of thermodynamic singularity partition func- 
tion is dominated by an exponential number of metastable states, while 



the ofF-equilibrium dynamics fails to reach equilibrium and falls into an 
asymptotic aging states with internal energy extensively higher then the 
equilibrium one [S]. The relaxation time to equilibrium is infinite at all 
temperature below the dynamical transition. As the Cortona meeting has 
testified, progresses have been recently achieved in the comprehension of 
both the thermodynamics and the dynamics of Mean-Field models. On one 
hand, the Guerra |H] and Talagrand "T, analysis shows that Parisi ansatz, 
which in these models take the simple one step form, describes well the 
thermodynamics, on the other Ben Arous, Dembo and Guionnet |S] could 
prove that the physicist's equations -and underlying assumptions leading 
to their derivations- can be fully mathematically justified. 

The situation is much less clear as one wants to deal with systems with 
finite range interactions. The low temperature properties of finite range 
spin-glass models, despite remarkable attempts of mathematical analysis 
in] remain an open issue of scientific discussion. While there is no consen- 
sus about the fate of mean-field thermodynamical spin glass phases when 
the range of interaction is finite, everybody would agree that the metasta- 
bility phenomena associated to the dynamical transition in ^spin models 
should become cross-overs as soon as the range of the interaction takes a 
finite range character. Possible scenarios for this cross-over have been put 
at the basis of a phenomenological theory of the glass transition known 
under the name of 'random first order transition scenario' which describes 
activated processes responsible for ergodicity restorations in terms of an 
effective droplet model where a bulk restoring force proportional to the 
configurational entropy competes with interface free-energy terms that op- 
pose to relaxation ^ E3 ^D ^] ■ Computations of the resulting effective 
barrier in the context of microscopic models have been presented in |13[I14| . 
Many researchers believe that an accomplished microscopic theory of this 
cross-over could provide important hints about the relaxation processes in 
supercooled hquids and the glass transition .15 . 

In order to understand the relation between mean field and finite range 
systems one can use the classical tools of Kac models JHIj where the in- 
teraction range R is considered a tunable parameter. One would like first 
to understand the equilibrium and dynamics in the Kac limit where the 
interaction range is sent to infinity after the thermodynamic limit |16l I17| , 
and then to study the dynamics for finite R and metastability effects in 
an asymptotic expansion in R. Recently a full mathematical control of the 
Kac limit for thermodynamics has allowed to prove convergence of free- 
energy and local correlations to the the mean-field values in p-spin models 
for even p |18[ [T^ . In this paper I would like to extend the analysis to the 
dynamics. In particular I will consider the case of the p-spin model and 
study: (a) the dynamics in the Kac limit (b) the form that the relaxation 
time should take in an asymptotic expansion in the range of interaction. 



This will lead to a picture where the long time relaxation can be seen as a 
passage from metastable state to metastable state, and the relaxation time 
computed as a decay rate. The analysis I will present will rely strongly on 
heuristic derivations. I believe that the analysis of the Kac limit, clear from 
the physical point of view, can be made fully rigorous with mild adaptation 
of the techniques of [5j. Conversely, the dynamics on R dependent time 
scales, remains to a large extent a challenge for theoretical physics, and the 
mathematical analysis seems to me far away in time. 

In several occasions in this paper I will refer to the dynamical and 
static properties of the mean-field p-spin model as they are known from 
the physicist's analysis. For many of these properties mathematical proofs 
are lacking. Throughout my reasoning I will ignore this fact and assume 
that the mentioned properties are correct. 

Some of the results presented in this paper have appeared in form in 
mi and US). 

2 The Kac Spherical j9-spin model 

In order to study the relaxation in a Langcvin setting, I use a variation 
of the locally spherical version of the Kac p-spin model considered in JS] • 
Consider A^ = {1,2, ■■■,L}'^ the hypercubic lattice of size L and periodic 
boundary conditions, partitioned into hypercubic boxes B^ of size £ {L/£ is 
a large integer), x G t^Lji- The model is defined in terms of real variables 
Si £ R subject to the local spherical constraint |^~'^X]ieB ^i ~ M <- ^ foi' 
some small epsilon and interacting through the random Hamiltonian 

(ii,...,ip)GAJ 

where the couplings Ji-^^...^i are for each p-uple of the lattice i.i.d. normal 
variables with zero mean and variance 

E{J!,....J = J2 -^M\^ - ^l\/R)■■■4'i\l - ^p\/R). (2) 

leAL 

The function ip : i?+ -^ i?+, is chosen to be positive ip{x) > and nor- 
malizable, /p ip{x)dx = 1. Notice that only variables that are at distances 
of order R or lower can effectively interact. The usual mean-field model 
is recovered choosing R = L and ip^x) = la;<i. I will consider instead the 
following regimes: (a) the Kac limit where L-^oo,R^oo,£-^oom the 
specified order, (b) the regime 

L» R»£» log(L) (3) 



where non trivial finite dimensional effects can be expected to take place 
|21|. Consider the Langevin dynamics 

where the ^^ (t) form a set of Lagrange multipliers chosen in such a way to 
enforce the spherical constraint in average at each time: E{£^'^J2i£B ^ii^)"^) — 



1, and rii{t) is white noise with amplitude \'2T and hi is an arbitrary extra 
field term that will be used to generate response functions. The equation 
will be interpreted in the Ito sense, and will be thought to be regularized 
at small times (e.g. discretized) so that a regular solution -differentiable 
with respect to h- will exist at all times. 

Below the temperature T^ = yj{p(jp — 2)p~'^ /(p — 1)^"^) the mean field 
model is in a non ergodic regime. One can expect for finite R the existence 
of two kind of time scales: time scales independent of i?, which should 
coincide with the ones of the corresponding mean field model and are inef- 
fective to relax at low temperature, and R dependent time scales on which 
there is ergodicity restoration. In order to understand the dynamics of the 
systems, it is interesting to consider its evolution in the Kac limit, where, 
in the order, L -^ oo, R -^ oo, ^ — > oo for fixed value of time. I would 
like to argue that for uniform initial conditions, the local correlation and 
response functions, defined respectively as 



ci"^''(i,") = ^E^^W'^^(") 



tend with probability one to non-random functions C(i, u) R(t, u) of the 
time arguments t and u, independent of x that verify the usual equations 
describing the dynamics of the mean field model (see e.g. :5^ for their 
detailed form and a review of mean-field off-equilibrium dynamics). 

Unfortunately, I cannot at the moment offer a proof of this statement. 
I would like however to discuss the indications that this is the case. The 
derivation parallels derivation in the usual mean field case, and shows that 
large R in the Kac case is analogous to large N in the mean field case. 

2.1 An exact equation 

In this section I derive an exact equation, valid for all values of the pa- 
rameters and all kind of initial conditions, that is a consequence of the 



Gaussian distribution of the couplings. Consider observables ^[S', ry] de- 
pending on the values of the variables and the Brownian noise at different 
times. Choices of interest will be A[S,r]] = Si{u) and A[S,t]] = i]i{u). I 
will denote with angular brackets (•) the average over the initial conditions 
and the thermal noise, and with E{-) the average over the quenched cou- 
plings Jij^^...^ip. From the Langevin equation, using partial integration on 
the Gaussian couplings one gets: 

E{^^A[S,v]) = -E{^,,{t)S,{t)A[S,v]) + {A[S,vMt)) + 

J2 EiAi,_,J'E—^ {A[S,r^]S^,,{t)...S,M (6) 

, ^^i.io .I'D 

In order to make explicit the J derivative, one can use the Martin-Siggia- 
Rose, Girsanov representation of the joint probability of paths Si{t) and 
Brownian noise r]i{t) starting from an initial condition S'i(O) chosen with 
probability fi{S{Oj) 






^ Ji,Z2,...,ipSi2{u)...Si^{u) + h^{u) + T]^{t) 
(i2,...Sp) 



±r,fit)]^iiSm (7) 



and observe, through integration by part, that insertion of iSi{u) in corre- 
lation functions of the kind {B[S, 77]), which improperly will be denoted as 
{iSi{u)B[S,ri]), acts as differentiation with respect to hi{u), or in turn as 
insertion of -^riiiu) 

{zS,{u)B[S,ri]) - ^A_(s[5,^]) = _L(s[5,^]^^(«)) (8) 



so that, 



Ejrj^ {A[S,liS.M■■■S^M - 

(-''Ji,i2,---Ap 



dvE^ — {r^,{v)S.,{v)...S.^{v)A\S,iT\S.,{t)...S,^{t)) 

+ Y,{S,{v)S,M■■■rl^Av)■■■S^MA[SAS^2{t)■■■S^M] 
1=2 J 



2.2 The ofF-equilibrium case 

If one choose an initial distribution independent of the quenched couphng, 
so that aj = 0, the last term in eq. lEl is zero. We assume in 

this section that this is the case, and that in addition /i is a translation 
invariant measure. One can consider for instance the uniform measure on 
the set I = {5* e i?^' : |£-i T,^(^B. Sf - 1| < e Vx}. 

Defining the empirical correlation function and response: 

i 

Ri{t>u) ^^^Y.^{\l-^\/R)S,{t)n,{u) (10) 

i 

and 



^^Y.^{\l-i\/R)pj\vEl^^{r^,{v)A[SMCi{t,vY'') 



1 

+ {p-l){S,{v)A[S,i^]Ri{t,v)Ci{t,vY-^)) (11) 

This is an exact equation valid for all values of L, _R, and t. From this 
equation one sees that if for typical realization of thermal noise and initial 
conditions, in the Kac limit the empirical functions Ci{t,u) and Ri{t,u) 
are self-averaging and tend to space homogeneous (^-independent) limiting 
functions C{t,u) and R{t,u), then the Lagrange multipliers would be x 
independent and one would have: 

E{S,{t)A[S,f^]) = -pi{t)E{S,{t)A[S,r^]) + E{A[S,r^Mt)) + 
^l dvE(J^{v,{v)A[S,r^])C{t,vr-' (12) 

+ {p-l){S,{v)A[S,7i])R{t,v)C{t,vr-^) . 

where I have used -^ J2i V'd^ — *l/^) ~^ 1- The being vahd for any A, it 
implies, in law, the effective, single site Langevin equation: 

S{t) = -fi{t)S{t) + rj{t) + e(t) + (p - 1) / dv C{t, v)P-^R{t, v)S{v) (13) 







where £,{t) is a Gaussian field with correlations {£,{t)(,{u)) = pC{t,uY ^, 
and subject to the consistency conditions: 



C{t,u) = {S{t)S{u)) 
1 



R{t,u)^—{S{t)^{u)). (14) 



where the bar indicates average over the distribution of the process. This 
process coincides with the one found in mean field dynamics and predicts 
that while at high temperature the system equilibrates after a finite time, 
at low temperatures T < Td the system fails to equilibrate and enters an 
off-equilibrium aging regime where the energy remains extensively higher 
then the equilibrium one. 

As soon as one gets away from the Kac limit and the range of interaction 
is finite, the relaxation time to equilibrium should always be finite. What 
the Kac limit is telling to us, is that for large enough R the relaxation 
time is essentially R independent at high temperature, while it becomes R 
dependent for temperatures T^T^- 

A proof of self-averaging and existence of the limit for the empirical 
correlation and response can be envisaged according to the lines of |S] 
for the mean field case. The main formal difference with respect to that 
case is that the white average over sites that is responsible for the self- 
averaging properties of the correlation and response functions in mean- 
field is here substituted by the weighed average ifTUIl . The self- averaging 
property, with existence of the limit, should appear when the number of 
variables participating to the average diverges. Space homogeneity follows 
from the hypothesis of statistical homogeneity of the initial condition. 

For the analysis that will follow it is interesting to investigate the behav- 
ior of the fluctuations for finite L, R and (.. of correlation and response on 
the scale i of the boxes Bx. Equation H13I) suggests that for fixed times i, s, 
typical fluctuations of correlations and responses © are of order 0{i~'^/^). 
Assuming the existence of a finite correlation length beyond which fluctua- 
tions become effectively independent, extreme value statistics implies that 
the maximum fluctuation is of order 0{^i'^/ \og{L)). The hierarchy ^ 
between length scales guarantees that with appropriate choice of L,R,£, 
this maximal fluctuation can be made arbitrarily small. 

2.3 Equilibrium dynamics 

To study the largest time scales in the system we need to understand equi- 
librium dynamics. As in the previous section, let us consider the Kac limit. 
The main difference with the off-equilibrium case is that the distribution 
of the initial state 

KSo) = ^e-^"^^o) (15) 

now depends on the quenched couplings. In the derivation to the analogous 
of equation ((T^ one should keep into account that 

^^^^- = S.,...S.^-(5,,...^.Je, (16) 



dJ, 



this gives rise to an additional term into the exact equation (|13|) for E{Si {t)A[S, rj]), 
that then reads: 

E{S,{t)A[S,v]) = ^^iAt)E{S^{t)A[S,v]) + E{A[S,iiMt)) + 

l^Y.iJ{\l-^\/R)pJ\vE(^■^{rJ.{v)A[S,rJ]Cl{t,vr-') + 

ip-l){S,{v)A[S,ii]Ri{t,v)Ci{t,vr-^)^ + 
■^ J2 H\l - ^\/R)E{S,{0)A[S, v]Q{t, Of-^) 

/,Z2,...,'ip I r—2 

E{{A[S,7j]S,,{t)...S,^{t)){SA,...S,^)eg) (17) 

In order to write closed equations in this case, besides hypothesizing the 
self-averaging condition for Ci and Ri discussed in the previous section, one 
needs to know the structure of the equilibrium correlations that appear in 
p7|l . The results of Q3 imply that the static correlations on scales of order 
R, tend to the corresponding mean field function in the limit. I am inter- 
ested here to the regime Tk < T < Td where these static correlations are 
vanishingly small in the limit, so that the last term in 1)17(1 will be zero. In 
equilibrium conditions the fluctuation dissipation theorem and time trans- 
lation invariance holds and R{t,s) = R{t — s) = (3 — q~^ ■ The resulting 
equation for C reduces after a little algebra to the usual equilibrium mean 
field equation |2"2j : 



TC{t) + ^ / dsC^-\t-s)—^ (18) 



dC{t) p [' 1 dCis) 

The behavior of this equation is very well known. Above Td the correlation 
function decays in a finite time to its equilibrium vale limt^oo C{t) — 0, 
while below Td, ergodicity is broken and limj^oo C{t) = qea- Qea is the 
Edwards- Anderson parameter, which does not coincide with the static value 
of the correlations in the Kac Hmit hmi<;^oo jd J2ieB^ Ei^i)lquii — ^■ 

We are in a situation where the Kac limit does not commute with the 
large time limit. A system with a finite but large R, for a large time will 
follow closely the evolution dictated by fl^ . but, eventually, will be able 
to relax below qea on R dependent time scales. Increasing R, one can 
make the relaxation time scale as long as wanted, in particular, it can be 
chosen so long to let the time to the system to explore ergodically the 
region where the local correlations Cx{t) never go below a certain value 
9o < Qea- This is a crucial observation that allows to identify typical 



equilibrium configurations as belonging to metastable states as it will be 
discussed in the next section. 



3 R dependent time scales 

On time scales diverging with the interaction range, the independence be- 
tween different spins which is implied by the self-averaging property of 
the correlation function should fail and the mean-field dynamical equa- 
tions loose their validity. In order to investigate the relaxation below qea 
one should proceed in a different manner, reminiscent of the Lebowitz and 
Penrose |23 analysis of metastability in Kac models of first order phase 
transition. 

I said that is a consequence of the fact that the ergodic time is divergent 
with R, for R large enough the equilibrium dynamics has the time to explore 
ergodically the region where for all boxes x the correlation Cx{t) > qo- This 
condition defines this region as metastable state in a sense that extends 
the one first specified in 'ST in connection with Kac models with first 
order phase transitions. In fact one has the following properties: 1) The 
system finds itself in states of constrained equilibrium. 2) The time to 
leave these states is large, but much larger is the time to come back once 
departed. 3) The correlation function is homogeneous in space. Formally, 
for typical equilibrium configurations S^, one can define the set TZ{Sa) = 
{S G T\ qx{S^S^) > q'^} and the metastable state as the constrained 
equilibrium state 

t^niSo){S) = ^^e-^^(^)l^(So)(5). (19) 

Notice that the states defined in this way do not form a partition of the 
equilibrium manifold. In fact couples of states so defined, corresponding 
to two different reference equilibrium configurations can overlap if these 
configurations are close enough to each other. Despite this fact, one can 
profitably describe the relaxation as a passage from metastable state to 
metastable state, and estimate the relaxation time as the inverse decay 
rate of these states. Denoting as H the Fokker-Plack operator and as 
TZ{So) the set complementary to TZ{So), the decay rate of TZ{So) is: 



A = / _ dSdS' {S'\n\S) X finiSo)iS)- (20) 

JSeTZ{So); S'€n{So) 

This expression can be bounded from above noticing that the Fokker- 
Planck operator is local, so that if S' has to be in TZ{So) then S has to be 
on the border of TZ{So), 

dUiSo) = {Se n{So)\ 3a;o : qxo{S,S'') = g"}. (21) 



Moreover, for any S, the integral / dS' (S"|7i|S') can be bound by a constant 
C and one gets the estimate 

X<C f dS M7^(s„) iS) = C^p^^ (22) 

where C is a constant of order 1. two comments are here in order: 1) 
Defined in this way, the escape rate A, and its estimate z^(S°) ^^^ func- 
tions of the reference configuration Sq. One can expect however that both 
quantities are self-averaging, i.e. assume values independent on Sq with 
probability approaching 1 in the thermodynamic limit. This suggests to 
consider 

Xtyp ~ exp [E{\ogZQniSo) - log^7j(5o))] = exp{~pAF{qo)) (23) 

where E stands here for Boltzmann-Gibbs average over the 5*0 and quenched 
average over the Gaussian couplings. 2) The definition of TZ{So) and of the 
decay rate depend on the value qo: which defines a sort of amplitude of the 
state. The "right" value optimizing the estimate, is the one such that a lo- 
cal equilibrium fluctuation of size go has equal probability to be reabsorbed 
in TZ{So) and to drive the system out of it. This, together with the fact 
that what we get is a lower bound for the decay rate, suggests that analo- 
gously to first order transition kinetics 23 , one should maximize AF{qo) 
with respect to qq. 

We get a relation between the relaxation time and a free-energy barrier 
defined coupling the system with a reference configuration. On general 
ground one can expect the barrier to be proportional to the interaction 
volume R'^. In the next section I will discuss some theoretical techniques 
to compute the barrier. 



4 The computation of the barrier 

I would like to discuss here possible strategies for the computation of the 
barrier. This is a hard task and my discussion will be highly conjectural. 
One needs to compute, ior S — TZ or S — dTZ, the following free-energy: 

Fs = ^{{J dSoe-^"^'o)^^gJ dSe'^His)^niSo)iS))) (24) 

where I denoted with ((•)) the average over the quenched couplings Ji^^...^i . 
This free-energy will consist in an extensive term plus corrections. In con- 
sidering the difference, the extensive term and the corrections due to finite 
size effects should compensate in the difference and only the term, finite in 
the thermodynamic limit, which is related to the relaxation time survives. 



10 



From the previous analysis, it is clear that F-ji is dominated by the config- 
urations where in all the boxes indexed by x one has that qx[S, So) = Qea- 
On the other hand, the restricted partition function Zg-jirgg\ can be written 
as a sum over all sites xq such that (72;o(S', Sq) — qo, namely, 

Xo 

Zx„ = f dSe-P"('^ (26) 

One can expect Z^g to be dominated by a single configuration consisting 
in a localized excitation of linear size 0{R) around the point xq. For 
large R, for most sites xq these excitations will become self-averaging and 
independent of the site xq- The partition function Zq^k^s^-^ will receive 
contributions from these sites, but could also receive contributions from 
exceptional sites particularly keen to excitation and where the relaxation 
would be initiated with the highest probability. At present I do not have 
a clear cut argument to decide if the dominant contribution to the barrier 
is given by typical or exceptional sites. In favor of the first hypothesis 
one could argue that the condition (, » logL should guarantee that too 
exceptional sites are absent from typical samples. If this is the case, any 
site Xq is equally likely to initiate the relaxation independently of Sq. In 
favor of the second one, one could argue that relaxation below qsA proceeds 
through amplification of small 0(l~'^l'^) fluctuations in the initial condition. 
This kind of amplification effect has been observed in realistic short range 
glass models, where sites that exhibit faster then average relaxation in the 
fast regime, are more likely to initiate the slow part of the relaxation |24| . 
In that cases however the interaction range cannot be considered a large 
parameter. 

At any rate, the free-energy Fq-jk^So) can be estimated through the the 
replica method, taking advantage of the identity 

E{FoniSo)) = lim -r((Z"-i I dSQe-f'"^'^\ZV,^^So)is,^ ~ !)/«»• (27) 

As usual, the computation is performed starting from integer values of m 
and n, and continuing to real values using a an appropriate modification 
of Parisi ansatz to take into account the constraints. It should be possible 
as in case of the computation of the unconstrained free-energy, to relate 
the replica computations to a variational principle in the space of "random 
overlap structures" |25| . but in this paper I will not pursue this path. 
Instead, I will sketch the outline of the replica computation of the free- 
energy in the simplest possible approximation |2(il 1181 1^ . For integer n 



11 



and m the average over the quenched noise can be explicitly performed. 
One has to introduce local overlaps between replicas 

Q"'''=£lE^^"^? «'^=1' ■■•'" + "^ (28) 

where the replicas from 1 to m—1 describe the normalization Z"^~^, replica 
number m describes Sq and replicas from m + 1 to m + n describe config- 
urations in TZ{So). In terms of these order parameters one can write the 
replicated partition function as 



Zrepi ^ lidx^o 1^ „„ Pge-^^M (29) 

where F[Q] is a coarse grained replica free-energy [TO] 

F[Q] = R'^SiQ] = R'' / cf^x [if ({Q„,/3}, x) + V{Q{x))] (30) 

Jo 

with 

K{Qo...p,x) = ^^[/(Q„,^(x))-/(g„,Mx))] 

a.J3 

V{Q) = -fE/W"^/3)-^Trlogg (31) 

and we have defined 

QaA^) ^ dyi){x~y)Qa^l3iy). (32) 



Notice that thanks to the scaling F[Q] = W^SIQ] for large R one can 
evaluate the integral H29(l by saddle point. I would like here just briefly 
discuss the results of the simplest ansatz for the matrix Qab{x), obtained 
assuming replica symmetry. It is clear from (|27|) that the integer n and m 
partition function is invariant under the group of S"^~^ x S'" of permutation 
of the m — \ replicas (a = 1, ..., m — 1 in our notation) coming from the 
denominator among themselves and the n {a = m + 1^ ...,n + rn) replicas in 
97?.(S'o) among themselves. Within the replica symmetric framework one 
chooses Xq = xq for all a > m. This ansatz for Xq is the translation into 
the replica formalism of the hypothesis discussed above that independently 
of Xq all sites are equally likely to start the relaxation. Replica symmetry 
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also implies that following structure of the matrix Qab{x) 

Qa,bix) = < 



1 


a~b 




six) 


a^b a,h — 1, ...,m— 1 




p{x) 


a < m, b > m 01 a > m, 


b < m 


q{x) 


a = m, b > m or a > m, 


b — m 


r{x) 


a ^ b a,b > m 





(33) 



It can be seen |2S] that the equation for s{x) decouples from the other 
variables and is solved by s{x) — p{x) — 0. As far as the other parameter 
are concerned, a further simplification consists in considering saddle points 
with r{x) = q{x), which can be verified to exist ^21- With this ansatz, the 
various terms of the free-energy become, in the limit of small n: 

K{{q{x)},x) ^ -n^[qP{x) ^ qPix)] 

V{q{x)) = -n (^^qP{x) + ^ [log(l - q{x)) + q{x)]^ . (34) 

and the action functional 

5= f dx[K{x) + V{x)] (35) 



V{q) is a single minimum function at high temperature (T > T^), while it 
has two minima below. The free-energy difference between the minima is 
well known to represent the configurational entropy per spin 5](r) in the 
mean-field model, i.e. the logarithm of the multiplicity of ergodic com- 
ponents divided by the volume |26|. Within the simple approximation in 
which the replica matrix is parametrized in term of a single space depen- 
dent parameter, we are in presence of Landau like field theory analogous 
to the one for systems with a first order transition. The free-energy barrier 
can the be estimated through instantonic techniques, looking at the mini- 
mum action solution with boundary condition equal to qEA at infinity and 
to go in xq mE|- Close to Tk this results in a conventional droplet theory 
where the free-energy barrier results from an effective droplet model where 
a bulk free-energy, in this case proportional to the configurational entropy, 
competes with a boundary term, which for the nature of the theory is just 
a surface term. The detailed computation of the barrier J3| leads to a 
Vogel-Fulcher like expression 

A^ = i?"^-^^ (36) 

with 7rf = d — 1 and C{T) a regular function of T. This computation 
supports in the context of a microscopic model the phenomenological "en- 
tropic droplet" arguments first introduced in ^ and recently revived in 
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|12|. Formula H36|l is obtained in the replica symmetric approximation and 
could be changed by considering better saddle points. About this result I 
have mixed feelings: on the one hand one finds a Adam-Gibbs like relation 
of inverse proportionality between configurational entropy and free-energy 
barrier, which is a good thing. On the other, in the relation one finds an 
exponent 7^ = d — 1 which would predict a lower critical dimension of 1 for 
the ideal glass transition and would be equal to 2 in three D. This is at vari- 
ance with the expectation, based on the Adam-Gibbs theory and verified 
in numerical simulations of model glasses, of a linear relation between the 
free-energy barrier and inverse configurational entropy in three dimensions. 
I believe that while the general form AF = R'^ -^Sd^ only depends on the 
fact that the bulk free-energy difference between the minima is the config- 
urational entropy, should be robust to improvements of the calculation, the 
value of 7(j could be modified by better approximations. One can envisage 
two sort of modifications in the replica computation: 1) Within the ansatz 
that Xq is independent of a, to find better saddle points for the matrix 
Qa,b{x). This has been attempted in ^3] where an replica broken solution 
to the instantonic problem has been found. Unfortunately, this solution 
though changing the function C{T), it does not change the exponent. It 
is my opinion that saddle points in this class cannot modify the exponent. 
2) One can consider more complex configurations where the position Xq of 
the center of the instanton depends on a. Saddle points of this kind would 
break the replica symmetry in a vectorial way |^, and could lead to a 
different value of the exponent 7^. The search of new saddle point to the 
replica variational problem will be matter of future research. 



5 Conclusions 

In this paper I presented the beginning of the analysis of the dynamics 
of the spherical p-spin model with Kac kind of interactions. I argued in 
favor of the existence, both in off-equilibrium and equilibrium dynamics at 
low temperature, of two different time regimes for large R. There is a first 
regime where dynamical time scales are independent (or quasi-independent) 
of R. Evidence for this regime is reached observing that in local correla- 
tions and responses should verify the dynamical mean-field equations of the 
infinite range model in the Kac limit. I am confident that the dynamical 
equations in this regime will be fully mathematically justified. Then there 
is a second dynamical regime, on time scales diverging with _R, where the 
system manages to explore ergodically the configuration space. The ex- 
treme separation of scales that one has for large R allows the description of 
the dynamics as a passage between metastable states that are explored er- 
godically before being abandoned. Thanks to that property, the relaxation 
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time can be estimated as the escape time from typical states and related 
to a well-defined free-energy barrier. The computation of this barrier is 
unfortunately a highly non-trivial problem. We presented here a simple 
attempt based on replica symmetry predicting an ideal glass transition for 
the mode in dimension d > 1 and giving a modified Adam-Gibbs relation 
between free-energy barrier and configurational entropy. The simple the- 
ory we have sketched could fail due to several kinds of replica symmetry 
breaking effects. While we believe that the general relation will continue 
to be valid, the value of the exponent in modified Adam-Gibbs relation 
and then the lower critical dimension for the ideal glass transition could be 
changed in the ultimate theory. 

Acknowledgments 

This work was supported in part by the European Community's Human 
Potential programme under contract "HPRN-CT-2002-00319 STIPCO". 

References 

[1] M. Mezard, G. Parisi and M.A. Virasoro, Spin glass theory and be- 
yond, World Scientific, Singapore (1987). 

[2] D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35 1792 (1975) 

[3] E. Gardner, Nucl. Phys. B297, [FS14] 74 (1985). 

[4] T.R. Kirkpatrick, P.G. Wolynes, Phys. Rev. B 36, 8552 (1987); T. R. 
Kirkpatrick, D. Thirumalai, P. G. Wolynes, Phys. Rev. A 40 (1989) 
1045 

[5] L. F. Cugliandolo in "Slow relaxations and nonequilibrium dynamics 
in condensed matter" Lecture notes, Les Houches, July 2002, Ed. by 
Jean-Louis Barrat et al. 

[6] F. Guerra Comm. Math. Phys. 233 (2003) 1 

[7] M. Talagrand, C. R. Acad. Sci. Paris, Ser. I 337 (2003) 111, Ann. 
Math. 163 221 (2006) 



[8] G. Ben Arous, A. Dcmbo, A. Guionnct, Preprint |m'ath.PR/0409273| 
To appear, Prob. Th. Rel. Fields (2006). 

[9] CM. Newman, D.L. Stein, J. Phys.: Condens. Matter 15, R1319- 
R1364 (2003) and references therein. 

[10] P. G. Wolynes, Jour. Res. NIST, 102 187 (1997), X. Xia, P. G. 
Wolynes, Proc. Nat. Acad. Sci. 97, 2990 (2000), X. Xia, P. G. Wolynes, 



15 



Phys. Rev. Lett 86, 5526 (2001), V. Lubchenko and P. G. Wolynes, J. 
Chem. Phys. 119 9088 (2003) 

[11] G. Parisi, "Gauge Theories, Spin Glasses and Real Glasses" Talk pre- 
sented at the Oskar Klein Centennial Symposium, cond-mat/9411115 
Proceedings of the XIV Sitges Conference "Complex Behaviour of 
Glassy Systems" 10-14 June 1996. Wolrd Scientific, Singapore 1996. 



[12; 

[13; 

[14 

[15; 

[16; 

[17; 
[is; 

[19 
[20 

[21 
[22; 

[23 

[24; 

[25; 

[26; 



J.-P. Bouchaud and G. Biroli, J Chem Phys. 121 7347 (2004) 

S. Franz, J. Stat. Mech. P04001 (2005) 

M. Dzero, J. Schmalian and P. G. Wolynes, preprint 
|cond-mat/0502011| 



For discussion of this point see e.g. M. Mezard, Physica A 306 25 
(2002) 

M. Kac, Phys. Fluids 2, 8 (1959) 

J.L. Lebowitz, O. Penrose, J. Math. Phys. 7 (1), 98-113 (1966). 

S. Franz and F. L. TonineUi Phys. Rev. Lett. 92 030602 (2004), J. 
Phys. A: Math. Gen. 37 (2004) 7433 

S. Franz and F. L. TonineUi J. Stat. Mech. P01008 (2005) 

S. Franz Europhys. Lett. 73 492 (2006) 

J.L. Lebowitz, O. Penrose, J. Stat. Phys. 3 211 (1971) 

A. Barrat "The p-spin spherical spin glass model" preprint 
Lcond-mat/9701031, unpublished. 

J.S. Langer, Ann. Phys. 41 (1967) 108. Ann. Phys. 54 (1969) 258. 

A. Widmer-Cooper and P. Harrowell, Preprint c ond-mat/0512035| D. 
Ashton and J.-P. Garrahan to appear. 

M. Aizenman, R. Sims, S. L. Starr, Phys. Rev. B 68 (2003) 214403 

S. Franz, G. Parisi, J. Phys. I (France) 5 (1995) 1401, Phys. Rev. 
Lett. 79 (1997) 2486, Physica A 261 (1998) 317 



[27] Vik. S. Dotsenko, M. Mezard, J. Phys. A 30 (1997) 3363, Vik S Dot- 
senko J. Phys. A: Math. Gen. 32 2949 (1999) 



16 



